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In this paper, we have studied Brans-Dicke Cosmology in anisotropic K ant owski- Sachs space- 
time model; considering variation of the velocity of light. We have addressed the flatness problem 
considering both the cases namely, (i) when the Brans-Dicke scale field cj> is constant (ii) when (j> 
varies, specially for radiation dominated era perturbatively and non-perturbatively and asymptotic 
behaviour have been studied. 

PACS numbers: 



I. INTRODUCTION 

In recent years theories with varying speed of light (VSL) has been attracted consid- 
erable attention [1-10] due to its ability to solve the so-called cosmological puzzles - the 
horizon, flatness and Lambda problems of big-bang cosmology. All previous attempts 
to overcome these difficulties invoke the basic inflationary form, where the observable 
universe experiences a period of 'superluminal' expansion. Here, one has to modify the 
matter content of the universe such that ordinary Einstein gravity becomes repulsive and 
halts the exponential expansion. One can however resolve to theories with varying speed 
of light as the alternative method to solve the above mentioned cosmological puzzles. 
Here, instead of changing the matter content of the universe, one has to change the 
speed of light in the early universe. The universe is assumed to be radiation dominated 
in the early stage and the matter content is assumed to be the same as in the standard 
big-bang (SBB) model. So the geometry and expansion factor of the universe go in 
accordance with the SBB model. Here for a free falling observer associated with the cosmic 
expansion, the local speed of light will decelerate from a very large value to its current value. 

So far some works have been done with VSL by Magueijo and co-workers [1-5] and others 
[6-10] specially in isotropic cosmological models (which also include BD-thcory). Very 
recently, Magueijo has investigated [11], the possibility of black holes formation by studying 
spherically symmetric solutions with VSL. In this paper, we will investigate in details the 
anisotropic Kantowski-Sachs (KS) model with VSL theory. The paper is organized as 
follows: section II deals with the basic equations in BD-thcory for KS model with varying 
speed of light. Perturbative solution to the flatness problem will be discussed in section III, 
while non-perturbative study of this problem will be discussed in section IV for radiation 
era only and the exact solution to the flatness problem is presented in section V. In section 
VI, we will study solutions to the quasi-flatness problem and examine its asymptotic 
behaviour. The lambda and quasi-lambda problems and their asymptotic behaviours will 
be investigated in section VII. The paper concludes with a short discussion in section VIII. 



II. THE BASIC EQUATIONS 

For anisotropic KS model with metric ansatz 
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ds 2 = -c 2 dt 2 + a 2 (t)dr 2 + b 2 (t)dft 2 , 
the BD field equations with varying speed of light are 
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and the wave equation is 
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Here the velocity of light c is an arbitrary function of time , u> is the BD coupling parameter 
and the BD scalar field is <j> — ^. From the above field equations we have the "non- 
conservation" equation 
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Now p = ^pc 2 is the equation of state in the radiation era for which the general solution 
is (from equation (3) ) 
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The above field equations in BD-theory have been formulated in Jordan frame. To switch 
over to Einstein frame we shall have to make the following transformations 
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dt=y/G4>dt, a^^/Opa, b=y/G$b, < 7 =( UJ +2) /n ( G( ^' P={G(p)~ 2 p, 

p = (G0)- 2 P , (6) 
and the above field equations become (c is treated as constant) 
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and 
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with the prime " ' " and " . " represent the differentiation with respect to t and t 
respectively. 

One can interpret these field equations as standard KS equations with constant G and a 
scalar field a is added to the normal matter . The scalar field behaves like a "stiff" perfect 
fluid with equation of state 
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If the velocity of light is constant , then in Einstein frame total stress-energy tensor is 
conserved but there is an exchange of energy between the scalar field and normal matter 
according to the following equation 
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On the other hand , if the velocity of light varies then we have two separate "non- 
conservation" equations 
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Thus standard KS model field equations with constant G for VSL are also the field equa- 
tions for BD theory in Einstein frame but here one must add to normal matter the scalar 
field energy and pressure . Therefore , the total energy density and pressure are given by 
p t = p + p a and p t =p + p a respectively . 



III. PERTURB ATIVE SOLUTIONS TO THE FLATNESS PROBLEM 



For a possible solution to the flatness problem in VSL theory , let us first study solutions 
when there are small deviations from flatness . The critical energy density p c in a BD 
universe is given by the equation 
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In Einstein frame the critical energy density for normal matter is 
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Hence, if we define the total critical energy as 



Pa= Pc+ Pa , 
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then from (15) we have 
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Let us define a relative flatness parameter as 
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whose evolution equation is given by 
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Here 7 — 1 = p/pc 2 = p/pc 2 , is the equation of state, with 7 constant. Since 
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in Jordan frame, the deviation from flatness is given by 
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where e = (p — p c )/pc measures natural deviations from flatness. Now 5, the adaptation 
parameter satisfies the differential equation 
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If 8 is assumed to be very small compare to unity i.e.,<5 <g; 1, then neglecting the square 
term in (21), one can integrate the above to get 
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where 5q is the integration constant. Hence we have 
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In BD theory G oc 4> 1 , so to solve the flatness problem cf) should decrease in the early 
universe and the first term on the right hand side must dominate over the second one. 

In particular, in radiation era (7 = 4/3) the above equation (23) simplifies to 
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where we have used equation (5) to obtain (f> . If a is small and c — cq [a™?) 2 '™^ 1 )] ™^ 3 ™ 2 -* 
then [12,3] it is still possible to solve the flatness problem for n < — 1 — y/7/3 . Furthermore, 
since p c oc l/(afe 2 ) 4 / 3 at late times, so the second term is always very small compare to the 
first and hence may be neglected. 
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IV. NON-PERTURBATIVE SOLUTIONS WITH <f> = CONSTANT: MATTER 
AND RADIATION DOMINATED ERA 

When <p = ^0 (constant) we have 5 = e in (20) and the differential equation for S, (21) 
can be written as 



e = £ (l + £ ) ( 7 - + 2( 7 -l)M +2-e (25) 
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Now, if we assume a power-law form in the scale factors of the velocity of light then the 
above equation becomes 
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In order to get an exact analytic solution to the above equation, we assume e <C 1 , so we 
can neglect e 2 -term. We then have the solution 

e = eoa 7+2n 2 /(3™-2) fo 2( 7 -l)+4n(«-l)/(3n-2) ^7) 

with e as the integration constant. 

Furthermore, without assuming any restriction on the parameter e, if we assume the two 
arbitrary constants n and 7 to be equal then also we have an exact analytic solution for e, 
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where A is an integration constant 
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V. EXACT SOLUTIONS TO THE FLATNESS PROBLEM: RADIATION 

DOMINATED ERA 

To obtain the exact solutions to the flatness problem we assume the variation of light 
mathematically by the relation [3,12] 
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However, in the Einstein frame the above relation reduces to c = c 
and equation (12) becomes 
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The equation of state p = (7 — l)/5c 2 simplifies the above equation further to 
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FIG. 1: Here the flatness parameter e in (27) has been plotted against t for the radiation era 
(7 = 4/3) choosing different values for the constant n. The scale factors here are assumed to be in 
simple expanding form. 
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FIG. 2: This figure shows the variation of the flatness parameter from (28), where there is no 
restriction on e. This figure is also drawn for the radiation era and assuming as before the simple 
expanding form for the scale factors. The variation of e for different choice of the arbitrary constant 
A is shown in the figure. 



An exact integral can be obtained (assuming n = 7) as 
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Hence for the radiation dominated era , we have in Jordan frame 
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where we have used p(ab 2 ) 4 ^ 3 = p(ab 2 ) 4 ^ 3 . 



Moreover, in the Jordan frame we have equation (29) for the velocity of light and we can 
write 
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Thus, if is a decreasing function then \c/c\ < \[na/a + 2(n— 1)6/6]/ (3n — 2)| and we can 
therefore conclude that a varying speed of light in the early universe (with weaker gravity) 
helps to solve the flatness problem. 

Alternatively for general equation of state p = (7 — l)pc 2 , if we make the transformation 



x = ^a 7 & 2(7_1) 

and assume the variation of the velocity of light as 
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then the equation of continuity (4) becomes 
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Integrating this equation gives 
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VI. SOLUTIONS TO THE QUASI-FLATNESS PROBLEM 

A. When BD scalar field 4> is constant 

In this section we shall investigate whether it is natural to have evolution which asymptotes 
to a state of expansion with a non-critical density . Let us start with the continuity equation 
for constant <f> and assume as before the equation of state 



P 

pc 2 



-7-1 



and the velocity of light to be 



c = c 
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then the integral gives 
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(It is to be noted that in order to obtain the integral one has to assume n = 7). 
Now substituting this value of p in equation (2) we have 



(38) 
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In particular, for radiation era the above equation simplifies to 
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If O is the density parameter, then it can be defined as 
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We note that the ratio between the two terms on the right hand side is almost a constant 
for a quasi- flat open universe. However, for 7 = 4/3, we get 

= 1 + 3^ (42) 
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For expanding universe a, b grows with time. So asymptotically the second term on the 
R.H.S will be negligible and we have 

o = -l 

3 

which leads to an open universe with finite O value today. 



B. When BD scalar field varies 

For the radiation era (7 = 4/3), let us define 

y = 0(ab 2 ) 2/ \ z = 0(a 2 b) 2 / 3 , c = c z m (43) 
So from the field equation (2) we get ( using equation (5) ) 
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Now, from equation (38) the expression of p takes the form 
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In addition the critical density p c is obtained from the differential equation 
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Hence the expression for the density parameter is 



q = P_ = Z V z + 3 y V 1 + 3 ) y 2 



n OV' z' z' 2 / i , 2ui\ a 2 



or equivalcntly, 
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Thus, if 2m + 1 < then as t — > oo, z — > oo so we have f2 — ► 1. But if 2m + 1 > then 
t — > oo, O — > —2m. So asymptotically, we have a quasi-flat open universe. 



C. General asymptotic behaviour 

In this section, we consider general asymptotic behaviour when 2m + 1 > i.e., for 
quasi-flat open universe. The equation (44) can be approximated by 
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For the case of fioo < 1, we have </> — > </>o, a(r) <~ r V n oo and 6(r) <~ t-V^oo i n the limit as 
t — ► oo , or equivalently, 

0-0o, a~r 1 /(«~+ 1 ), 6^ T V(n»+i) ) as t ^ oo. 

In addition, when £1^ = 1, we get a ~ t 1 / 2 , 6 ~ i 1 / 2 , as expected for flat radiation 
asymptote . 



VII. THE LAMBDA AND THE QUASI-LAMBDA PROBLEMS 

In this section, we shall examine the effect of variation of speed of light when a cosmological 
term is introduced in the BD field equations. In fact, the incorporation of a cosmological 
term is equivalent to introduction of a vacuum stress which obeys an equation of state 
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Thus equation of continuity can be generalized to 
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Also the field equation (2) can now be written as 



(49) 



Now, using the above expression for vacuum stress-energy and using 7 — 1 = p/ pc 2 as the 
equation of state for the matter and assuming 
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the above equation of continuity can be integrated (assuming n = 7) to give 
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On substituting this value of p in the above field equation (49), we have 
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If we define the generalized density parameter as 
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If we now substitute the value of p from equation (50) we obtain 
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when 7 > 1. The 2nd, 3rd and 4th terms are negligible compare to the first term, hence 
we have 
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So, for large a, b (i.e., a,b — > oo) 



Furthermore, if A = then 
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So for the radiation era (7 = 4/3) 
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As a, 6 — > 00, fi/(f2 — 1) — > 8/14 i.e., O = —4/3. This asymptotic behaviour is same as in 
Sec. VI A where for quasi- flatness problem the BD scalar field is assumed to be constant. 



VIII. DISCUSSION 



In this work an extensive analysis of the BD solutions for anisotropic cosmological 
model has been done where there is a variation of the velocity of light with time we 
have assumed the velocity of light to be in power-law form in the sacle factors. The 
flatness problem has been discussed in details. Here perturbative, non-perturbative and 
exact solutions of flatness problem have been obtained. The graphical representation for 
non-perturbative solution in the radiation era has some interesting feature. Quasi-flatness 
problem with general asymptotic bahaviour has also been discussed and in most cases 
we obtain an open universe. Finally, we have also studied the lambda and quasi-lambda 
problem in VSL and the asymptotic behaviour is very similar to that in the flatness problem. 
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